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ABSTRACT 

We consider the merging of compact binaries consisting of a high-mass 
black-hole and a neutron star. From stellar evolutionary calculations which 
include mass loss we estimate that a ZAMS mass of ^ 80 Mq is necessary before 
a high-mass black hole can result from a massive 0-star progenitor. 

We first consider how Cyg X-1 with its measured orbital radius of ~ 17 Rq 
might evolve. Although this radius is substantially less than the initial distance 
of two 0-stars, it is still so large that the resulting compact objects will merge 
only if an eccentricity close to unity results from a high kick velocity of the 
neutron star in the final supernova explosion. We estimate the probability of 
the necessary eccentricity to be ~ 1%; i.e., 99% of the time the explosion of a 
Cyg X-1 type object will end as a binary of compact stars which will not merge 
in a Hubble time (unless the orbit is tightened in common envelope evolution 
which we discuss later). 

Although we predict ~ 7 massive binaries of Cyg X-1 type, we argue 
that only Cyg X-1 is narrow enough to be observed and that only it has an 
appreciable chance of merging in a Hubble time. This gives us a merging rate of 
~ 3 X 10~^yr~^ in the galaxy, the order of magnitude of the merging rate found 
by computer driven population syntheses, if extrapolated to our mass limit of 
80 Mq ZAMS mass for high-mass black hole formation. Furthermore, in both 
our calculation and in those of population syntheses, almost all of the mergings 
involve an eccentricity close to unity in the final explosion of the 0-star. ^From 
this first part of our development we obtain only a negligible contribution to our 
final results for mergers, and it turns out to be irrelevant for our final results. 

In our main development, instead of relying on observed binaries, we consider 
the general evolution of binaries of massive stars. The critical stage is when the 
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more massive star A has become a black hole and the less massive star B is 
a giant, reaching out to A. We then have a common envelope, and we expect 
hypercritical accretion to A. A will accept a small fraction of the mass of the 
envelope of B but will plunge deep into B while expelling B's envelope. 

We expect that star B can at least be in the mass range 15 ~ 35 Mq while 
the black hole A has a mass of 10 Mq. About 20 percent of the binaries of this 
type are found to end up in a range of orbital radii favorable for merging; i.e., 
outside of the relevant Roche Lobes, but close enough so that these final binaries 
of compact objects will merge in a Hubble time. The narrow black-hole, 0-star 
orbits do not seem to be found in population syntheses because in them mergers 
happen almost completely as a result of kick velocities. In the exception. Case 
H of Portegies Zwart & Yungelson (1998) which includes hypercritical accretion, 
common envelope evolution is more effective and we are in agreement with their 
results. 

We find that the high-mass black-hole, neutron-star systems contribute 
substantially to the predicted observational frequency of gravitational waves. 

We discuss how our high mass for high-mass black hole formation can 
be reconciled with the requirements of nucleosynthesis and indicate that a 
bimodal distribution of masses of black holes in single stars can account, at 
least qualitatively, for the many transient sources which contain high-mass black 
holes. 
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1. INTRODUCTION 

The supernova community used to believe that stars above a certain mass, about 
30 — 40 Mq ZAMS (zero age main sequence) will collapse into a massive black hole (MBH) 
of mass of order 10 Mq. The argument was that in these stars the mantle was bound with 
a binding energy well above 10^^ erg so that the supernova shock was not strong enough to 
expel it. 

Whereas this may be true for single stars, Woosley, Langer & Weaver (WLW,1995) 
showed that in binaries, where the hydrogen envelope of the primary star has been 
transferred to the companion in RLOF (Roche Lobe Overflow), the evolution of the 
resulting "naked" He star; i.e., the star without hydrogen envelope, led to a substantially 
smaller presupernova core than that of a single star with hydrogen envelope. A comparison 
of compact core masses from naked He stars and those evolved by Woosley & Weaver (1995) 
for single stars is shown in Figure 1, taken from Brown, Weingartner & Wijers (1996). 
Detailed reasons for the great difference in the evolution of "clothed" and naked He cores 
are given in WLW (1995). 

Stars with ZAMS masses > 40 Mq lose their masses by strong winds, whether in 
binaries or not, and become Wolf-Rayet stars. In an earlier paper WLW (1993) investigated 
ZAMS masses of 35, 40, 60 and 85 Mq. In those up through 60 M© the hydrogen envelope 
was blown off early enough for the He cores to evolve as naked ones and compact core 
masses were around 1.5 Mq (gravitational). In fact, with inclusion of extensive mass loss 
the lower line in Fig. [l|, which heads just above 1.5 Mq for the higher ZAMS masses, 
gave the WLW (1993) correspondence of Fe core mass to ZAMS mass for single stars of 
masses 35 — 60 Mq. Thus, rapid mass loss by wind in this mass region which removes the 
H-envelope before appreciable He core burning begins leaves a He core which burns as a 
"naked" one. (In the case of the massive stars the relation shown in Fig. |l| between Fe core 
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mass and He core mass no longer holds because of large wind losses in the latter.) In the 
case of the 85 Mq star some hydrogen envelope remained during an appreciable part of the 
He burning, so the He core burned as a (partially) clothed one and the compact core was 
more massive, in the range 1.7 — 2.0 M©, depending upon ^^C(q;,7)^^0 burning rate. For 
the very massive stars, the He core burns, at least part of the time, as if clothed. The mass 
loss rate used by WLW (1993) has been shown to be too high. Such rates were obtained 
from the free-free fluxes and modelling of infrared spectral lines which involve the density 
quadratically. Measurements from the polarization of the radiation (St-Louis et al. 1993; 
Moffat & Carmelle 1994) obtain a mass loss rate of < 50% of that used by WLW. The 
polarization depends linearly on the density, provided that the wind is optically thin to 
electron scattering. Thus, estimates based on this are independent of the inhomogeneities 
which are known to be present. The lower mass loss rates come into agreement with the 
M from dynamical arguments. Certainly the WLW calculations should be redone with the 
lower rates, which may change our conclusions. 



2. THE COMPACT STAR 

Thorsson, Prakash & Lattimer (1994) and Brown & Bethe (1994) have studied the 
compact core after the collapse of a supernova, assuming reasonable interactions between 
hadrons. Initially, the core consists of neutrons, protrons and electrons and a few neutrinos. 
It has been called a proto-neutron star. It is stabilized against gravity by the pressure of 
the Fermi gases of nucleons and leptons, provided its mass is less than a limiting mass Mpc 
(proto-compact) of ~ 1.8 Mq. 

If the assembled core mass is greater than Mpc there is no stability and no bounce; the 
core collapses immediately into a black hole. It is reasonable to take the core mass to be 
equal to the mass of the Fe core in the pre-supernova, and we shall make this assumption. 
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although small corrections for fallback in the later supernova explosion can be made as in 
Brown, Weingartner & Wijers (1996). If the center collapses into a black hole, the outer 
part of the star has no support and will also collapse. We then get a massive black hole 
containing the entire mass of the presupernova star, perhaps of order <^ 10 Mq. (At least in 
binaries the companion star will have removed the hydrogen envelope by either Roche Lobe 
Overflow or by common envelope evolution, depending on the mass of the companion, and 
wind will diminish the He core before the SN explosion.) 

If the mass of the core is less than Mpc, the electrons will be captured by protons 

p + n + u (2.1) 

and the neutrinos will diffuse out of the core. This process takes of order of 10 seconds, 
as has been shown by the duration of the neutrino signal from SN1987A. The result is a 
neutron star, with a small concentration of protons and electrons. The Fermi pressures of 
the core are chiefly from the nucleons, with small correction from the electron. On the 
other hand the nucleon energy is increased by the symmetry energy; i.e., by the fact that 
we now have nearly pure neutrons instead of an approximately equal number of neutrons 
and protons. Throsson et al. (1994) have calculated that the maximum mass of the neutron 
star Mns is still about 1.8 M©; i.e., the symmetry energy compensates the loss of the Fermi 
energy of the leptons. Corrections for thermal pressure are small (Prakash et al. 1997). 

The important fact is that the ten seconds of neutrino diffusion from the core give 
ample time for the development of a shock which expels most of the mass of the progenitor 
star. 

But this is not the end of the story. The neutrons can convert into protons plus K~ 
mesons. 



n ^ p + K 



(2.2) 
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Since the density at the center of the neutron star is very high, the energy of the K~ is very 
low, as confirmed by Li, Lee & Brown (1997) using experimental data. By this conversion 
the nucleons can again become half neutrons and half protons, thereby saving the symmetry 
energy needed for pure neutron matter. The K~ , which are bosons, will condense, saving 
the kinetic energy of the electrons they replace. The reaction (|2.2|) will be slow, since it is 
preceded by 

e- ^ K~ + u (2.3) 

which is actually effected in the star by reactions such as e~ + p ^ + p + u, with 
strangeness breaking. (Times are long enough for chemical equilibrium to be realized.) It 
becomes energetically advantageous to replace the fermionic electrons by the bosonic i^~'s 
at higher densities. Initially the neutrino states in the neutron star are filled up to the 
neutrino chemical potential with trapped neutrinos, and it takes some seconds for them 
leave the star. These must leave before new neutrinos can be formed from the process ( |2.3| ). 
Thorsson et al. (1994) have calculated that the maximum mass of a star in which reaction 



( p.2|) has gone to completion is 

Mnp ^ 1.5 Mq , (2.4) 

where the lower suffix NP denotes their nearly equal content of neutrons and protons, 
although we continue to use the usual name "neutron star". This is the maximum mass 
of neutron stars, which is to be compared with the masses determined in binaries. The 
masses of 19 neutron stars determined in radio pulsars (Thorsett & Chakrabarty, 1998) are 
consistent with this maximum mass. 

The core mass Mc formed by the collapse of supernova must therefore be compared to 
the two limiting masses, Mpc and Mj^p. If 

(/) Mc>Mpc (2.5) 
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we get a high mass black hole, of mass essentially equal to the full mass of the presupernova 
star. If 

(77) Mpc > Mc> Mnp (2.6) 

we get a low mass black hole, of mass Mq- Only if 

(777) Mc < Mnp (2.7) 

do we get a neutron (more precisely, "nucleon") star from the SN. Only in this case can we 
observe a pulsar. In cases (II) and (III) we can see a supernova display. In case (I) only 
initial neutrinos from electrons captured in the collapse before Mq becomes greater than 
Mpc can be observed but no hght would reach us. 

We tentatively choose the lower limit of ZAMS mass for making MBHs to be 80 Mq. 
On the other hand it is beheved that ZAMS above 100 Mq do not exist, because of excessive 
formation of electron pairs. So we assume that the range of 80 to 100 Mq is available. 



3. RATE OF FORMATION 

We are interested in massive binaries containing one star of ZAMS mass between 80 
and 100 Mq. As in Bethe & Brown (1998) we start from the assumption that there is one 
supernova per century per galaxy in a binary. Assuming also that 10 Mq is required for a 
star to end up as a SN of type II (or lb or c), the formation of binaries of M > 10 Mq is 
also 10~^ per year per galaxy. We assume a Salpeter function with index n = 1.5; then the 
fraction of such stars between 80 and 100 Mq is 

/ = 8-^/2 - 10^^/2 = 1.26 X 10"^ . (3.1) 

So the rate of formation is 

ai = 1.26 X 10~V~Ver galaxy . (3.2) 
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We require that this star A be accompanied by a companion S of M > 10 Mq. 
Assuming that the distribution of mass of the companion is dq, with 

q = Mb/Ma (3.3) 

and having assumed Ma — 90 Mq on average, we need q > 1/9 which has a probabihty 

1 - g - 0.9 , (3.4) 

hence formation rate 

q;2 = (1 - q)ai = 1.13 X lO^V"^ ■ (3.5) 

In order to observe strong X-rays from the MBH, the distance a between star A, the 
MBH, and star B, an O or B star, should not be too large, let us say < 150 R©. On the 
other hand it must not be < 30 R© because otherwise the two stars would merge already 
at this stage of evolution. Assuming, as in Bethe & Brown (1998) a distribution da /I a, the 
probability that a falls in the desired limits, is 

p = 7"Mn(150/30) = 0.23 (3.6) 

giving for the probability of formation (per galaxy) 

as = 0.23q;2 = 2.6 x 10"V~^ ■ (3-7) 

Although we take the same logarithmic distribution as in Bethe & Brown, we think of 
the lower limit as a > 2 x 10^ km so that the more massive stars we deal with here he 
inside their Roche Lobes and the upper limit as a < 2 x 10^° km, although the latter is 
uncertain because 0-star binaries probably will not be recognized for such a large separation 
(Garmany et al. 1980). High mass transfers from the O-star occur only when it nearly fills 
its Roche Lobe. It is then bright for a time (Massevich et al. , 1979) 

T = 2.7 X 10^ yr . (3.8) 
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so the expected number of strong X-ray emitting binaries in the galaxy is 

AT = agr = 7 . (3.9) 

Only one such binary has been observed, Cygnus X-1, at a distance of 2.5 kpc. It 
might be thought that X-ray binaries at larger distance have escaped detection because of 
absorption of X-rays in the highly ionized galactic medium. But in fact, X-rays of energy 
> 3keV should be seen throughout the galaxy, penetrating even the Galactic disc. On 
the other hand, we are probably seeing the closest binary of Cyg X-1 type, with orbital 
separation only about half of the initial separation in the double 0-star progenitor. Thus, 
other such objects may well have substantially lower luminosities. The 0-star in such a 
binary would probably be in the disc and might well not be observed at greater distances, 
so that the X-rays could not be associated with the binary. In any case, it is clear that in 
the disc there is only one bright Cyg X-1 type object, a fact we use below. 

In the Magellanic clouds with less than one-tenth the mass of the galaxy, two high-mass 
black- hole binaries, LMC X-3 with a B-star companion and LMC X-1 where the donor is 
probably an 0-star, have been observed. We believe the high incidence of HMBHs in the 
LMC to be not only a consequence of the large amount of star formation there, but also of 
the low metallicity compared with solar. We return to a discussion of high-mass black hole 
formation in Section |8|. 



4. MERGERS 

In the case of low-mass black holes (LBH), the rate of mergers is limited chiefly by the 
disruption of the binary due to the recoil that star B experiences when it goes supernova. 
For massive black holes (MBH) this is of little concern since the orbital velocity of the 
binary before the SN process is already of order 600 km s~^ (see below), comparable to the 
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higher recoil velocities, so that the recoil is unlikely to distrupt the binary. 

Instead of this, the concern for MBH is the two stars of the binary may be too far 
apart so that they will not merge during Hubble time (assumed to be 10^° years). Note that 
the progenitor binary must have separation a > 30Rq in order that the two 0-stars not fill 
their Roche Lobes. At present, the separation in Cyg X-1 is only slightly more than half 
this, which is permissible because now only the 0-star needs to be in its Roche Lobe, the 
black hole has negligible radius. But as we show below, even this small separation is still 
large enough to substantially cut down the probability of merging within a Hubble time. 
Because of rapid mass loss in such massive stars of ZAMS 80 — 100 Mq there is no well 
established procedure for calculating their evolution in binaries. Starting from a separation 
> 30 Mq mass exchange by the very massive star with the lower mass companion when 
the former is on the main sequence (Case A mass transfer) could decrease the separation 
substantially, perhaps causing the stars to coalesce. However, the resulting binary would 
then widen (roughly back to its original separation) as the Wolf-Rayet, which remains after 
mass transfer from the very massive star, loses mass by wind. We do not have any evidence 
of the net change (S. Portegies Zwart, private communication). 

Because of the uncertainties in the evolutionary scenario, we adopt an empirical 
approach, beginning from the measured separation a, in Cyg X-1. The requirement that 
the two stars merge in a Hubble time then sets a lower limit on the eccentricity e of the 
final neutron-star black- hole orbit: 

a2(l-e^)<ao (4.1) 

where 02 is the semi-major axis of the orbit after the SN event and e is the eccentricity of 
the orbit; qq will be calculated in Section ^ 

The left hand side of ( |4.1| ) is related to the angular momentum J of the post-SN orbit. 



In fact 

J = a^Vsil - e'y/' (4.2) 

where V3 is the orbital velocity if the orbit were circular. (We do not need to consider the 
mass of the system.) By Newton's laws, 

aalf = G{Ma + Mn) (4.3) 

where Ma is the mass of the MBH and M„ that of the neutron star resulting from the SN 
of star B. Hence 

.f = G{Ma + M„)a2(l - e^) . (4.4) 

On the other hand, immediately after the SN, the distance between stars A and B is still 
ai, and the relative velocity is 

V2 = Vi + Q (4.5) 

— * _ 

where Q is the recoil velocity due to the SN. The angular momentum is 

J = ai X V2 

J2 = a^lfsin^^ (4.6) 

where ip is the angle between V2 and ai. 

This angle can be calculated but is complicated. We now make the approximation of 
replacing sin^ ip by 

(sin2^) = ^. (4.7) 



Then, equating ( |4.4| ) and (^]6|) and using ( |4.1| ) 



al Vi < ^G{Ma + M„) ao . (4. 
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Eq. ( [4 .71) involves a double approximation; (1) is uniformly distributed over the sphere, 
and (2) only the average of siv? ip is used, instead of calculating the limit on aiV2 separately 
for each possible direction of Q and then averaging the results. We believe that a correct 
calculation would give a somewhat higher limit on aiV2- 

The product aiVg^ can be calculated assuming that the orbit of the binary is circular 
before the SN event, an assumption which is probably very nearly correct. Then, similar to 

o, 

aiV^ = G{Ma + Mb) (4.9) 

where Mb is the mass of star B, the companion of the MBH, before the SN event. In 
Cyg X-1, it is believed that the mass of the MBH is = 10 M0 and that of the 0-star 
companion M'^ = 17 M© (Herrero et al. 1995). After RLOF of the 0-star these masses 
happen to be just exchanged, assuming conservative mass transfer, so that immediately 
before the SN, 

Ma = 17 Mq, Mb = 10 . (4.10) 

However, we shall use "standard" values 

Ma = Mb = 10 M0. (4.11) 

The SN changes B to a neutron star, M„ = 1.4 M©, while A remains unchanged. We have 
now 

al Vl = G{Ma + Mb) ai. (4.12) 

^Later we shall show that common envelope evolution of the black hole in the evolving 
0-star is more likely. 
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Dividing (gj) by (CT) 

^2 3 / Ma + M„ \ ao 



Usually, ai > so V2 < Vi. /^From ([4 .51) this means that the recoil velocity Q of the SN 
must be approximately opposite to the pre-SN orbital velocity Vi, and of the same order of 
magnitude. Setting 

Q-V^i = -QVfx (4.14) 
Q/Vi = to (4.15) 



thus (|4.13| ) become 



1 +uj^ -2ujfi< B (4.16) 

r-.<'-^. (...) 

Since B is in general considerably smaller than 1, we thus get the condition 

stating once more that Q must be of the same order as Vi. The fraction of the total solid 
angle available to Q is (1 — /i)/2. 



Now according to ( [4.12| ) 



, 2 10~^ X 20 X 2 X 10^3 1012 

= - X = 2.7 X 10^^ X (4.19) 

3 tti ai 

Vi = 520km s"^(10^V«i)^^^- (4-20) 

In Cygnus X-1. the distance between the two stars is 

ai = 17Rq = 1.2 X lO^^cm. (4.21) 

So ai is likely to be of order 10^^. This means Vi is of order 500 km s~^, and Q is 
likewise. Now Cordes & Chernoff (1997) have found that the recoil velocities of pulsars are 
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distributed bimodially with cr = 175 and 700 km s ^, respectively. The smaller velocity 
Gaussian is negligible for our Q. For the larger cr, we set 

X = Q/a = 1. (4.22) 

The Gaussian distribution is 

f dx = 27r-^/2^2^xe-^'. (4.23) 

We assume Vi = a for the Gaussian with the greater a, then for this Gaussian, u = x and 
the probability is, cf. ( [4.17|) 

/ -^[5 - (1 - xy]f dx^'^ B'/^ = 0.14 B'/\ (4.24) 
J ix o Y vr 

For large B, substitute 0.95 for B^^'^. Considering that only 20% of pulsar recoils are in 
the larger a Gaussian, the fraction of binaries containing massive black holes which merge 
during a Hubble time is 

= 0.03 (4.25) 
with B given by ( |4.13| ). Inserting our values for M^, Mb and M„, 

B = 0.85 ao/ai. (4.26) 

It remains to calculate Cq. 



5. MERGERS, NUMERICAL 

If the orbit is circular, merger will take place within a Hubble time if the initial distance 
is < Rmax, with Rmax giveu by Bethe & Brown (1998), eq. (7.7) 

= 1-6 X 10^^ km^ -. (5.1) 
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R is the same as the semi- major axis introduced in Section |, so Oq = Rmax- If the orbit is 
excentric, the merger time for a given major axis is diminished by a factor (P. Eggleton, 
1998, private communication) 

Z(e) = (1 - e2)3-689-0.243e-0.058e2 _ (5 2) 

Therefore the condition on the initial semi-major axis, 02, is 

4 Z{e) < C. (5.3) 

Taking the fourth root 

a2 Z^''' < R^a. (5.4) 
Z^/^ is almost 1 — e^, so the left hand side of (5^) is almost 02(1 — 6^), thus 



02(1 - e^) < i?„,,(l - e2)0.0777+0.0608e+0.0195e2_ (5 5) 

The last factor depends only slightly on the eccentricity. So we replace in it by an 
average of 0.5. This makes that factor 



Inserting in (|5]1|) Ma = 10 Mg, Mb = 1.4 Mq, 

i?max = 7.1x lO^km. (5.7) 



The upper limit of 02(1 — e^) is the definition of Oq, according to ( |4.1| ). Thus the right hand 
side of ( |5.5| ) equals oq , and we have 



ao = 6.5 X 10^ km (5.8) 

and (|4.13| ) becomes 

5 = 5.1 X 10^ km/ai. (5.9) 
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The measured distance between the two stars in Cygnus X-1 is 1.2 x 10'^ km. If this distance 
does not change in RLOF, then ai = 1.2 x 10'' km, 

B = 0.42 (5.10) 
= 0.009. (5.11) 

This estimate can be checked from Fig. 4a of Kalogera (1996). In our situation her 
parameters are (3 ~ 0.6, C, ~ 0.9. Using the Eggleton formula eq. ( p.2|) , we find that 
e > 0.76 for merger. The average survival probability is ~ 0.3 for these eccentricities and 
our ^ is ~ 0.3, so that the final probability is 1/5 x 0.24 x 0.3 ^ 0.014 close to our 0. 

Multiplying this by the formation rate as, eq.( p.7|) , we find for the rate of mergers 

= 2 X lO^V"^- (5.12) 

So this is a rare event as we might expect. In fact, (|5.12|) is probably still an overestimate 
because we predicted 7 X-ray emitting binaries in the galaxy (eq.( p.9|) ) and we see only 
Cyg X-1. In our assumed distribution of these objects out to 150Rq, the others would be 
expected to have substantially larger ai's than Cyg X-1, in which case their merger rate 
would be negligible. Thus, a possible guess of the rate of mergers is not much larger than 

1/7 of (Km. 



We compare our result with results from population synthesis: 

1. First we note that the smallness of the rate of mergers results from the separation at 
the black-hole, 0-star stage being too large for merger (for spherical orbits) without a 
high eccentricity from kick velocities. Whereas Portegies Zwart & Yungelson (1998) 
obtain a (high- mass) black- hole, neutron-star birthrate of 18.7 x 10^^ yr~^ in the 
galaxy without kick velocities, their merger rate is zero for this case. Lipunov et al. 
(1997) find an order of magnitude greater merger rate with inclusion of kick velocities 
than without. 
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2. Portegies Zwart & Yungelson use 40 Mq for their ZAMS mass for black hole limit 
and we compare with those for 40 Mq from Lipunov et al.. Increasing the mass limit 
to 80 Mq and taking into account the factor 1 — q would decrease the results of these 
authors by a factor of ~ 5. Thus, Portegies Zwart & Yungelson would have a merging 
rate of 2 x 10^^ yr^^ and Lipunov et al., one of ~ 4 x 10^^ yr^^ for our mass limit of 
80 Mq. Our rate obtained by taking 1/7 of ( |5.12D is 



R = 3x 10"^ yr"^ (5.13) 

about the same as the latter. We believe that our simple considerations give same 
understanding of the results obtained in the latter population synthesis, and we shall 
return to a discussion of the work of Portegies Zwart & Yungelson later. 

We can also try an "observational estimate" of the merging of Cyg X-1 type objects, 
although we have only one such object near its Roche Lobe. Assuming as in ( |3.8| ) Cyg X-1 
to be bright for a time r = 2.7 x 10^ yr, the birth rate of such bright objects 

= 3.7 X 10~^ jTs'\ (5.14) 

Given the merging probability of the resulting black-hole, neutron-star binary = 0.009 
from (|5.11|) we arrive at a merging rate of 4 x 10~^ yr~^ for the galaxy, in rough agreement 



with (|5.13|) . This may be an underestimate since Cyg X-1 will probably go into stable mass 



transfer after losing substantial mass to the black hole. LMC X-3 with roughly equal black 
hole and B-star masses (Kuiper et al. 1997), may already be in stable mass transfer. Thus 
the 2.7 X 10^ yrs may well be an underestimate. 

Although the last two sections give results for mergers which are negligible compared 
to our final results in §7, once hypercritical accretion is included in our common envelope 



^We compare with their Case H in which hypercritical accretion was included in the next 
section. 
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evolution, we believe that they are instructive in that they match rather well the results of 
population synthesis, especially those of Lipunov et al. In Portegies Zwart & Yungelson 
(1998) (generally unstable) mass transfer starts on the dynamical time scale of the donor, 
changing along the way to the thermal time scale and the last part of the mass is sometimes 
even transferred on a nuclear time scales (Portegies Zwart, private communication). Their 
initial separations of the 0-star binaries are nearly double ours. (We take ours before 
substantial main-sequence burning which increases the stellar radius by a factor ~ 2.) 
Thus, their common envelope evolution does tighten the orbits, but clearly not enough 
for substantial merger of the compact objects evolved later. The exception is Case H 
of Portegies Zwart & Yungelson, the case similar to ours with inclusion of hypercritical 
accretion, which gives a merger rate ~ 35 times that of Cases B-G. We develop in the next 
section why hypercritical accretion is so effective. 



6. COMMON ENVELOPE EVOLUTION 

In the last two sections we have addressed the evolution of massive binaries in which 
the primary evolves into a high-mass black hole in an empirical way. We match rather well 
the population synthesis results of Lipunov et al. (1997). In particular, these authors find 
that introduction of kick velocities increases their merging rate by an order of magnitude. 
Portegies Zwart & Yungelson (1998) find zero mergers without kick velocities ! We beheve 
this to result from the fact that their black-hole 0-star binaries lie too far apart for the 
resulting binary of compact objects to merge in a Hubble time without the increase in 
merger distance given by an eccentricity close to unity; i.e., these authors have very little 
(no) merging from circular black-hole neutron-star orbits, except for their Case H which 
includes hypercritical accretion which we discuss below. 

Clearly the population syntheses of Lipunov et al. (1997) do not include efficient 
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common envelope evolution which bring the compact object (high-mass black hole) close to 
the He core of the companion star. 

We now develop our scenario for common envelope evolution, which hypercritical 
accretion makes particularly plausible, in that we have quantitative control over the crucial 
quantities. As the envelope of the giant expands to meet the compact object one can see 
from Rasio & Livio (1996) that the common envelope evolution begins rather quickly. Soon 
after the beginning of mass transfer the compact object creates a tidal bulge in the evolving 
companion, transferring angular momentum to the companion. As the giant companion 
loses mass, the isentropic envelope responds by expanding. The compact object plunges 
into the companion, the chief loss in orbital energy occuring in ^ 1 year in time (Terman, 
Taam & Hernquist 1995; Rasio & Livio 1996). Our analytical estimate is roughly consistent 
with this because our hypercritical accretion proceeds at ~ 10^ Eddington, or about 1 Mq 
yr~^ and the total accretion is of the order of solar masses, giving the time in years. 

In the literature common envelope evolution is thought to occur when the two stars 
involved do not differ too much in mass, by not more than a factor of 2 to 3. This condition 
arises because, as higher mass ratios, the mass receiving star will then not be able to accept 
the large amount of mass transferred to it on the short thermal time scale of the companion. 
This thermal time scale is (van den Heuvel 1994): 

= 10^)/(M/Me)^) (6.1) 

In fact, the factor 2 to 3 in g initially came from a factor ~ 10 in rth, from the work of 
Kippenhahn & Meyer (1977) who considered Case A mass transfer (during main sequence 
evolution). In Case B mass transfer the transfer takes place as the giant traverses the 
Herzsprung gap, much more quickly than in Case A, so it is not clear why q could not be 
substantially reduced and still have common envelope evolution. 

The situation with hypercritical accretion is different from either Case A or Case B 
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(RLOF) mass transfer. Initially our ~ 10 Mq black hole is met by the expanding red giant 
or supergiant envelope. The black hole accretes some of the matter and transfers enough 
energy to the remaining matter to expel it. The convective envelope has constant entropy 
and must expand in order to replace the accreted and expelled matter. The black hole 
drops in gravitational potential closer to the He core of the companion in order to furnish 
the necessary energy to expel most of the matter. This whole process happens very rapidly. 
(In our case of hypercritical accretion, an accretion disc will be set up around the black 
hole, which cannot immediately accept the matter because of the high angular momentum 
of the latter. We assume the vicosity to be high enough so that angular momentum, but 
little mass, will be advected outwards. Observed masses ~ 7 Mq of high mass black holes 
in transient sources as discussed by Brown, Lee & Bethe (1999) substantiate this scenario.) 

Hypercritical accretion sets in when the envelope density reaches ~ 10^^ cm~^. From 
eqs. ( |5.3| ) and (|5.7| ) in Bethe & Brown (1998) we can show that the contribution to the 
coefficient of dynamical friction Cd from matter inside the accretion radius Rac-, i-e. from 
matter which can accrete onto the compact object, is 5cd = 2. Now the total coefficient of 
dynamical friction is q = 6 — 8 (Ruffert 1994; Ruffert & Arnett 1994). The remainder of q 
comes chiefly from the wake, at greater distances than Rac- Thus for 

Cd > 2 (6.2) 

we encounter highly nonconservative mass transfer and it is plausible that common envelope 
evolution ensues for a wider range of q than usually thought. In fact, as we show below, out 
of the ~ 20 Mq hydrogen envelope of a 30 Mq companion star only ~ 3 Mq is accepted by 
the black hole. 

[FTom the above argumentation we believe that q ^ 2 strongly favors common 
envelope evolution. (For Q = 2 the mass transfer is nearly conservative.) The companion 
star should have mass less than 35 — 40 Mq, because for higher masses it loses mass in 
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an LBL phase. It should have somewhat higher mass than the black hole, say 15 Mq in 
our case, so that mass transfer is unstable throughout the transfer of envelope mass but 
the final transfer is so rapid that it will overshoot into the range of stable mass transfer. 
As noted, Cyg X-1 is probably included in our interval of masses favorable for common 
envelope evolution. 

Having established the plausibility of common envelope evolution for the black hole 
in the expanding H envelope of an 0-star with ZAMS mass of 15 — 35 Mq, we now carry 
out this common envelope evolution for a typical 0-star mass of 20 Mq following Bethe & 
Brown (1998). The He core is 6 M©. 

We choose the coefficient of dynamical friction to be 6. In the Bethe & Brown 
notation Ma is the compact object mass, Mb the companion mass, and 

Y = Mea-K (6.3) 

The initial and final Y are related by 

.y^x 1+1/(0.-1) _ 2.4M^^ 

\yJ Ma, ^ ' ' 

depending only on the initial masses. The ratio of final to initial black-hole masses is given 

by 

so the final black-hole mass is 13 M© for an initial 10 Mq. We find 

ai _ Mb, Yf 



af MBf Yi 

where Mb^ is the companion He core mass of 5.3 Mq. 



12 (6.6) 
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7. MERGERS AFTER COMMON ENVELOPE 

The most important result of the calculation at the end of Section ^ is that ai/aj is 
large, of the order of 12. In ( |6.6|) , there are two factors; Yf/Yi is fairly large, about 4. 
Mb J Mb J is also fairly large, it is the ratio of the entire mass of the star to its He core 
taken to be (0.3)"^ by Bethe & Brown (1998). 

Thus the radius of the orbit shrinks by a large factor in the common envelope. We 
know the range of the final radius a/: its maximum is given by the condition that the 
merger should occur within a Hubble time. According to ( p.7|) , it would be about 7 x 10^ 
km if the orbit were circular. Eccentricity may raise this about 30%, 

^rnax _ g ^ 10"cm. (7.1) 

The minimum a/ is some multiple of the radius of the He star which star B becomes 
after removal of its hydrogen envelope. The He star radius may be about 5 x 10^°, and we 
estimate 

a™" = 1.5 - 3 X lO^^cm. (7.2) 

So aj^"^/a™" = 3 ~ 6. We assume the distribution of a to be da/7a, so the probability of 
having a/ in the permitted range is 

p = i^~^ = 0.15^0.25. (7.3) 

The center of the useful range is at about 5 x 10^^ cm. Multiplying by ai/af = 12, we 

get 

ai < lO^^'cm. (7.4) 

This is too large a distance to permit appreciable X-rays to come from the black hole. If 
the two compact stars are to merge by gravitational waves within a Hubble time, their 
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precursors (a black hole and an O- or B-star) can presumably not be observed as emitting 
X-rays. 

Conversely, objects like Cygnus X-1 probably will not lead to an observable 
gravitational- wave merger. When the 0-star expands as a giant and goes into a common 
envelope with the black hole, the orbit will contract to the extent that the black hole falls 
into the He core of the giant; there is a merger but its gravitational waves will be too low 
frequency to be observed. Thus the discussion in our Sections H and ^ is actually irrelevant, 
given our scenario for common envelope evolution. 

Returning to our discussion up to ( [7.3|) , we have to examine the values of the initial 
mass of star B. It should be greater than 15 Mq to permit common envelope evolution but 
less than 35 Mq so as to avoid undue mass loss by wind. With Maq = 90 Mq, this means a 
range of q of 0.2. Multiplying this by ai, in ( ^.2| ) and by the mean of p in ( |7.3| ), the rate of 
merger of a massive black hole with another compact star is 

^4 = 4 — 6 X 10~^yr~^ per galaxy. (7.5) 

This is in remarkably good agreement with the results for Case H with hypercritical 
accretion of Portegies Zwart & Yungelson, if we decrease their number by a factor of 5 to 
take into account our greater high-mass black-hole mass limit; namely, with this decrease 
they would have R = 7 x 10^^ yr~^ per galaxy. 

The great uncertainty is in ai because we do not know very well the minimum ZAMS 
mass which leads to a massive black hole, thus ai is likely to be uncertain by a factor 2 
either way. Together with p in (|7.3| ), the uncertainty of merger is about a factor 3, so 

04 = 1 ~ 10 X 10~^ yr~^ per galaxy. (7.6) 

This rate is much smaller than that of mergers of small black holes with neutron stars, 
as ~ 10~^, as one might expect from the difficulty of forming massive black holes. But 
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the signal to noise ratio in the gravitational wave detector depends on the chirp mass as 
(Mc/iirp)^'^- The chirp mass is 

M^mrp = /i°■'M°■^ (7.7) 

where /z is the effective mass and M is the total system mass. We have Mchirp ~ 3.3 Mq 
for our high- mass black- hole, neutron-star binaries, as compared with ~ 1.6 Mq for the 
low-mass black-hole, neutron-star binaries. Thus, effectively the high-mass black-hole 
merging rate should be multiplied by a factor of ~ 6. The two types can be distinguished 
by the chirp mass. 

Star B, after its supernova event, will be either a neutron star or a small black hole of 
mass slightly over 1.5 Mq, it has no occasion to accrete extra mass. Thus, we estimate that 
signals from massive black hole mergers with neutron stars should have a frequency about 
1/3 of those from smaller black holes. Taking the central value of 5 x 10^^ from ( [7.5[ ) and 
multiplying it by 6, we obtain an effective increase of ~ 30% over the merging rate of 10~^ 
yr~^ for the low-mass black-hole, neutron-star binaries. 



8. OTHER EFFECTS OF HYDROGEN "CLOTHING" 

We have seen that the He core of stars evolves quite differently according as it is 
"clothed" with an H envelope or not. We believe that this difference may explain problems 
in some recent investigations. 



8.1. Formation of black holes with light companions 

Portegies Zwart, Verbunt & Ergma (1997) have discussed the formation of black holes 
in low-mass X-ray binaries. Many of these transient X-ray sources have been discussed 
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recently with black- hole masses in the probable range ~ 6 — 7 Mq. Assuming a lower limit 
of ~ 40 Mq for black-hole formation, Portegies Zwart et al. find a much too low rate of 
formation. 

In their evolution, the black hole originates from the more massive component (A) of 
the binary. The less massive component (B) will spiral in to A when A becomes a giant 
and touches B. They show that B will survive this spiral- in only if its original distance 
from A is at least several hundred solar radii. Otherwise, B will spiral into the He core of 
A, will merge with A and thus be lost. In order for A to have a radius of several hundred 
Rq, it must have completed He core burning. This means that A has burned He while 
"clothed" with most of its H envelope. Consequently, the upper curve in Fig. 1 applies; the 
massive star A essentially burns as a single star. In this case stars of ZAMS mass as low as 
20 — 25 Mq can go into a black hole. We suggest this as a possible solution to the Portegies 
Zwart et al. underproduction of HMBH's. 

Our scenario also suggests an explanation of why the transient black holes are generally 
accompanied by low-mass companions. As the massive progenitor A of the black hole 
evolves as a giant, RLOF will generally transfer its H envelope to a massive companion 
(if one exists) during hydrogen shell burning or early in the He core burning phase. After 
RLOF the He core of the primary will burn as a "naked" He core, with possible fate as 
either low-mass black hole or neutron star (Brown, Weingartner & Wijers, 1996) but not as 
a high-mass black hole. Therefore a high-mass black hole will generally not have a massive 
companion, except in the relatively rare cases. 

If, on the other hand, star B has small mass, it can accept only very little mass from A 
in RLOF. Instead, interaction of the two stars will wait until A becomes a supergiant and 
its surface reaches B. Then, in the common envelope, B spirals in as described by Portegies 
Zwart et al. The He core of star A, mass about 1/3 of the original ZAMS mass, evolves as 
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described by Portegies Zwart et al., finally collapses as a SN and leaves a remnant of mass 
approximately equal to the He core mass, 6 — 7 Mq. 

8.2. Cygnus X-3 

Cyg X-3 may be a progenitor of a binary of two low-mass black holes. We believe the 
envelope of the 0-star progenitor of the compact object (A) in Cyg X-3 was probably lifted 
off in RLOF. This star A then had a SN event making it into a neutron star. When star 
B, the present He star, become a giant, star A would spiral in which would have converted 
the neutron star into a low- mass black hole (if it wasn't one already). The He star B in 
Cyg X-3 is certainly burning as "naked" . Our estimated is that the ZAMS mass of the 
progenitor of this He star is about the maximum mass that does not go into an LBL stage 
(because the envelope must have been used up in common envelope evolution), say 35 Mq. 
^Prom our fig. 1 the most likely fate of such a star, with H envelope lifted off in RLOF, is a 
LMBH, but a neutron star of mass ~ 1.5 Mq cannot be excluded. 

8.3. Nucleosynthesis 

The galactic ratio of oxygen to iron depends on the ZAMS mass above which single 
stars evolve into high mass black holes (cut-off mass) and therefore do not return matter to 
the galaxy. Oxygen is chiefly produced in quiescent burning before the supernova explosion 
and the amount is roughly proportional to ZAMS mass. Fe, on the other hand, is produced 
explosively in the SN explosion, and the amount is roughly independent of ZAMS mass, 
possibly decreasing slightly as the latter increases. Tsujimoto et al. (1997) have recently 
used the [O /Fe] ratio and observations in metal poor stars to determine a cutoff mass of 



Mcutoff = 50 ± 10 M, 



(8.1) 
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They remark that the influence of the metallicity dependence of stellar wind losses may 
be significant, but do not take it into account. We believe, however, the larger wind losses of 
stars with solar metallicity to have important effects. Over a wide range of ZAMS masses, 
between ~ 35 — 40 M© and ~ 80 M© according to Woosley, Langer & Weaver (1993), the 
winds remove the H envelope sufficiently rapidly that the He cores evolve as "naked" . Thus, 
the cores evolve into either low-mass black holes or neutron stars, their (gravitational) core 
masses following the lower line heading slightly above 1.5 Mq in Fig. 1. Brown, Weingartner 
& Wijers found that the primary in 1700-37 probably evolved into a low-mass black hole 
(I.e., it doesn't pulse.) whereas 1223-62 is known to contain a neutron star. In both cases 
the ZAMS mass of the primary was found to be ~ 40 Mq, so this approximately locates 
the ZAMS mass corresponding to M^rp of Section ^ in the case of solar metallicity. (As 
noted earlier, we estimate the ZAMS mass corresponding to Mpc to be ~ 80 M©.) It may 
well be that for metal poor stars, with much weaker winds, the hydrogen envelope is not 
removed rapidly enough for their He cores to evolve as "naked" . In this case the limit (|8.1| ) 
for ZAMS masses giving core mass Mpc could be more appropriate. 

We have not attempted to carry out a calculation of nucleosynthesis in our scenario, 
but wish to point out features of our scenario which will tend to increase the Mcutofi above 
which element production for the galaxy ceases. 

Disregarding mass loss, single stars above a certain mass Mmin evolve into high-mass 
black holes. From the upper curve in Fig. one might take the minimum He mass to be 
6 — 8 Mq corresponding to a ZAMS mass 20 — 25 Mq. But as the mass of the single star 
increases, there is big mass loss before the star ever reaches the supernova stage. So with 
extensive mass loss WLW(1993) find that ZAMS mass 35 Mq leads, after SN, to a compact 
object of about 1.5 Mq, thus a low-mass black hole. Thus high-mass black holes are only 
formed by single stars in a limited mass range, from 20 — 25 Mq to about 35 — 40 Mq. 
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Stars in this intermediate mass range, of width 10 — 20 Mq, do not return their matter, 
especially their Fe, to the galaxy. This group of stars is in addition to the stars above 
Mcutoff- To compensate for this fact, Mcutoff will have to be raised above the Tsujumoto 
value. Because the abundance of stars decreases with increasing mass, the raise has to be 
more than 10 — 20 Mq. Thus we suggest that Mcutoff may be as high as the 80 Mq which 
we used for other reasons in section |^. 

We note from the WLW (1993) calculations that "naked" He cores evolve into less 
massive carbon/oxygen cores than "clothed" ones. Thus we see that inclusion of mass loss 
will tend to move the mass above which all nucleosynthesis ceases even higher. Therefore, 
we believe that inclusion of the mentioned effects may decrease the apparent discrepancy 
between Meutofr of Eq. and our ZAMS mass of 80 Mq for making MBH's. 

We return to the higher observed incidence of HMBH's (LMC X-1 & LMC X-3) in the 
LMC. Suppose that the limit of Eq. ( |8.1j ) applies for the lower metallicity of the LMC; i.e., 
stars with ZAMS masses ~ 40 — 80 Mq in the LMC would not experience as large wind 
loss as those in the disc, because of the lower metallicity in the LMC. This could help to 
explain why two HMXBs containing high-mass black holes are observed in the LMC, and 
only Cyg X-1 is seen in the Galactic disc. 



9. CONCLUSIONS 

Our chief result is that those massive binary systems in which the primary is sufficiently 
massive > 80 Mq to go into a high-mass black hole can contribute importantly to observable 
gravitational waves upon merger of the final binary of compact objects. These mergers 
come chiefly from companion 0-stars in the range of 15 — 35 Mq ZAMS mass. 

We argue that the situation created by hypercritical accretion is favorable to common 
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envelope evolution in this range of companion masses, and that the time scale of this 
evolution is very fast. 

We show that population syntheses do not have a common envelope with efficiency 
comparable with ours, possibly because they generally do not include hypercritical accretion. 
Certainly they do not have as effective a tightening of orbits as we have. The exception is 
Case H of Portegies Zwart & Yungelson (1998) which does include hypercritical accretion. 
For the same mass limit for black hole, their results are in good agreement with ours. At 
first sight it may seem surprising that we match so well the Portegies Zwart & Yungelson 
work, provided they include hypercritical accretion as we do. Their calculation includes 
various phases of mass transfer before forming the binary consisting of a neutron star and 
black hole. However, results depend only on the ratio of the logarithmic interval favorable 
for gravitational merger with avoidance of coalescence to the total logarithmic interval over 
which binaries are distributed. Whereas the favorable logarithmic interval is shifted around 
in the various mass transfers, its magnitude is unchanged because the only scale is the 
radial separation of the two objects. 

Our final merger rate of ~ 4 — 6 x 10~^ yr~^ increases the Bethe & Brown (1998) rate 
of gravitational waves by a factor of ~ 1.3, largely because of the higher chirp mass with 
high-mass black holes. 

This important effect comes in spite of much higher mass limit ZAMS 80 M0 for 
evolution into high-mass black holes in binaries. We justify this high mass limit from results 
for the evolution of "naked" He stars by Woosley, Langer & Weaver. Namely, the more 
massive primary in binaries has its mass lifted off early, either RLOF or in LBL stage for 
the stars with ZAMS masses > 35 — 40 Mq. The resulting "naked" Wolf-Rayet or He star 
is deprived of the H-envelope which normally "insulates" He cores during their burning so 
the convective carbon burning stage is not skipped as it is in single stars of ZAMS mass 
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^ 20 Mq. The great entropy loss during the long duration of this stage results in a low-mass 
compact object. 

We show that the high mass limit for evolution into high-mass black holes in binaries 
is consistent with nucleosynthesis, because single stars with ZAMS masses in the range of 
20 to 35 — 40 Mq evolve into high-mass black holes without return of matter to the galaxy. 

We also indicate that we could evolve enough transient sources, binaries of a high-mass 
black hole and low mass main sequence star. 

There is interest in the merging of binaries composed of two high-mass black holes 
(Brady, Creighton & Thorne 1998). With our mass limit of 80 Mq for high-mass black 
holes the initial separation of the two massive 0-star progenitors must be > AORq for both 
of them to he inside their Roche Lobes. Case A mass transfer between two nearly equal 
mass stars can only widen the binary, since aj/ai = {ni/ufY and jj, will decrease in the 
transfer. Stars of ZAMS mass ^ 40 M0 lose their H envelopes rapidly in the LBL phase, 
so a common envelope phase in which the orbit is tightened is excluded. The Wolf-Rayets 
which result after loss of the H will lose most of their mass by wind, down to ~ 10 — 15 Mq 
before going into black holes. With spherically symmetric mass loss, af/ai = Mi/Mf ~ 4. 
Thus, the separation of the two final black holes will be several times the initial > 40 Rq 
which precludes merging in a Hubble time. We are in agreement with Portegies Zwart & 
Yungelson (1998) on this point. 

We beheve that in Brown & Bethe (1998) together with this paper we have given a 
consistent description of the evolution of most binaries of compact stars containing a black 
hole. 
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FIGURE CAPTION 

FIG. 1- Comparison of the compact core masses resulting from the evolution of single stars 
(filled symbols), Case B of Woosley & Weaver (1995), and naked hehum stars (Woosley, 
Langer & Weaver, 1995) with masses equal to the corresponding He core mass of single 
stars. The horizontal dashed lines indicate the mass of the heaviest known well-measured 
pulsar and the maximum mass of a neutron star. 
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Fig. 1. — Comparison of the compact core masses resulting from the evolution of single stars 
(filled symbols), Case B of Woosley & Weaver (1995), and naked helium stars (Woosley, 
Langer & Weaver, 1995) with masses equal to the corresponding He core mass of single 
stars. The horizontal dashed lines indicate the mass of the heaviest known well-measured 
pulsar and the maximum mass of a neutron star. 



